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Abstract 

Background: MicroRNAs are small RNA molecules that negatively regulate gene expression by translational 
inhibition or mRNA cleavage. The discovery that abnormal expression of particular miRNAs contributes to human 
disease, including cancer, has spurred growing interest in analysing expression profiles of these molecules. 
Quantitative polymerase chain reaction is frequently used for quantification of miRNA expression due to its 
sensitivity and specificity. To minimize experimental error in this system an appropriate endogenous control gene 
must be chosen. An ideal endogenous control gene should be expressed at a constant level across all samples and 
its expression stability should be unaffected by the experimental procedure. 

Results: The expression and validation of candidate control genes (4.5S RNA(H) A, Y1, 4.5S RNA(H) B, snoRNA, U87 
and U6) was examined in 21 rat cell lines to establish the most suitable endogenous control for miRNA analysis in a 
rat model of cancer. The stability of these genes was analysed using geNorm and NormFinder algorithms. U87 and 
snoRNA were identified as the most stable control genes, while Y1 was least stable. 

Conclusion: This study identified the control gene that is most suitable for normalizing the miRNA expression data 
in rat. That reference gene will be useful when miRNAs expression are analyzed in order to find new miRNA 
markers for endometrial cancer in rat. 

Keywords: Endogenous control genes, microRNA, Endometrial cancer 



Background 

Endometrial cancer is the most frequently diagnosed gy- 
naecological malignancy in the western world. According 
to the World Cancer Research Fund (WCRF) around 
288 000 women worldwide develop endometrial cancer 
annually and 74 000 died from this cancer in 2008 [1]. 
This malignancy can be divided into two different types, 
endometrioid adenocarcinoma (type I) and serous car- 
cinoma (type II). Type I develops from endometrial 
hyperplasia and is characterized by oestrogen depend- 
ence, well differentiated cells, and good prognosis. Type 
II, on the other hand, develops from atrophic endomet- 
rium and is oestrogen independent, with poorly dif- 
ferentiated cells, and has poor prognosis. The most 
common type, endometrioid adenocarcinoma (EAC), 
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accounts for approximately 75% of the reported cases 
[2,3]. As all cancers, EAC is a complex disease, where 
interactions between genes and environment play 
important roles in the course of the disease. Studies of 
such diseases are complicated due to the genetic hetero- 
geneity and environmental diversity of the human 
population, and therefore inbred animal model systems 
are often used to minimize the genetic diversity and to 
control the environmental factors [4] . The use of animal 
models does, however, not replace research on human 
material; these complement each other. 

Virgin females of the BDII/Han inbred rat strain 
spontaneously develop EAC at a frequency of more 
than 90% during their lifetime [5,6]. Due to the high 
degree of conservation between the human and rat ge- 
nomes, genetic aberrations in the BDII rat can be used 
to predict alterations associated with EAC development 
in humans. The similarities in physiology and patho- 
genesis between rats and humans imply that the BDII 
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rat strain is a suitable model for analysis of human 
endometrial cancer [7]. 

MicroRNAs (miRNAs) are small RNA molecules that 
regulate gene expression by inhibition of protein transla- 
tion or by degradation of their target mRNAs. It is esti- 
mated that up to 30% of human genes are regulated by 
miRNAs [8], and it has been shown that the functions of 
the target genes include a broad range of important bio- 
logical processes, such as development, differentiation, 
growth and metabolism [9]. As for most cancers certain 
miRNAs are differentially expressed in EAC compared 
to normal tissue (see ref. [10] for a review). However, to 
our knowledge, no such miRNA expression analysis has 
yet been done in models of EAC in rat. 

Due to its sensitivity and specificity quantitative poly- 
merase chain reaction (qPCR) is commonly used for 
miRNA expression analysis [11]. The data from qPCR can 
be analysed and presented as absolute or relative values. 
In absolute quantification a standard curve is used to 
calculate the quantity of the unknown sample, whereas 
relative quantification is based on the relative expression 
of the gene of interest compared to one or more reference 
genes [12]. For relative quantification methods, raw qPCR 
data are normalized in order to correct for variation 
caused by the amount of starting material, variation in re- 
action efficiency and sample purity. This normalization 
can be done using one or more stably expressed endogen- 
ous control genes. A gene that is used as a reference must 
have stable expression independent of sample and tissue 
type [13]. In most qPCR studies of miRNA expression 
small nuclear and/or nucleolar RNAs have been used as 
control genes, but there are also studies in which stably 
expressed miRNAs have been used. There are several 
reasons why small non-coding RNA genes might be well 
suited as controls for normalizing miRNA expression data, 
including their stable expression, such as the nucleotide 
size, and the use of identical assay chemistry [14,15]. 
Several human and mouse small RNAs have been tested 
and confirmed as suitable endogenous controls for quanti- 
fication of miRNA expression levels [14], but as far as we 
know no such study has yet been performed for small 
RNAs in rat. 

In this study, six endogenous control genes were 
selected and tested for expression stability in tissues 
derived from rat endometrium; five small nuclear RNAs 
and one small cytoplasmic RNA. The identification of 
suitable endogenous control genes is an important initial 
step in expression analysis since usage of an unstable 
control gene for normalization could result in mislead- 
ing conclusions. 

Results 

Differences between the triplicates were tested for signifi- 
cance by ANOVA (p < 0.05), which revealed no significant 



differences between the three replicates in any of the six 
endogenous controls. Students £-test were applied for 
differences in gene expression between malignant and 
non-/pre-malignant samples were tested for in each of the 
six candidate genes using the (Table 1). No significant 
difference in expression between malignant and non-/pre- 
malignant samples was detected for 4.5S RNA (H) A, 4.5S 
RNA (H) B, snoRNA and U87, but significant differences 
were found for U6 and Yl (p < 0.05). 

In order to identify and rank the most suitable control 
genes the data were analysed by two different algo- 
rithms, geNorm and NormFinder. GeNorm calculates 
the stability values (M) for all candidate genes, then 
eliminates the candidate gene with the highest M and 
repeats this procedure until the two most stable genes 
are left. The analysis revealed that all genes showed 
M-values below the geNorm default limit of 1.5 (Table 1). 
The control genes with the lowest M-values and highest 
expression stability were U6 and 4,5 S RNA (H) A, both 
with M = 0.714. U87 was the third most stably expressed 
gene with M = 0.797. NormFinder calculates a stability 
value (standard deviation, SD) for each gene, where a low 
SD represents more stable expression (Table 1). In con- 
trast to geNorm, Normfinder identified U87 and snoRNA 
as the two most stable genes with SD = 0.611 and 0.662, 
respectively. Although the two rankings differ, both algo- 
rithms identified Yl as the most unstable gene. 

When more than one control gene is identified by 
Normfinder, this algorithm will also calculate the accumu- 
lated standard deviation. The optimal number of control 
genes is indicated by the lowest value for the accumulated 
SD. For this data the lowest value is reached when five 
control genes are used and all of them are nuclear RNAs, 
whereas the excluded gene (Yl) is cytoplasmic (Table 2, 
Figure 1). 

Discussion 

It has been clearly demonstrated in numerous studies 
that miRNAs play critical roles in the development of 

Table 1 Statistical and stability analyses of candidate 
endogenous control genes by f-test, geNorm and 
NormFinder 



f-test* geNorm NormFinder 



Control gene 


P value 


Rank 


M 


Rank 


SD 


Rank 


Total i 


U87 


0.158 


3 


0.797 


3 


0.611 


1 


1 


snoRNA 


0.937 


1 


0.948 


4 


0.662 


2 


1 


4,5S RNA(H) A 


0.194 


2 


0.714 


2 


0.886 


5 


3 


U6 


0.023 


6 


0.714 


1 


0.696 


3 


4 


4,5S RNA(H) B 


0.373 


4 


1.016 


5 


0.773 


4 


5 


Y1 


0.032 


5 


1.152 


6 


1.255 


6 


6 



*The f-test refers to differences in gene expression between malignant and 
nonVpre-malignant samples. 
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Table 2 Panel of cell lines that were used in the present 



study 


Tumour 


Genetic background 


Pathology 


NUT6 


(BDIIxBN)xBDII 


EAC 


NUT43 


(BDIIxBN)xBDII 


EAC 


NUT50 


(BDIIxBN)xBDII 


EAC 


NUT81 


(BDIIxBN)xBDII 


EAC 


NUT128 


(BDIIxBN)xBDII 


EAC 


NUT48 


(BDIIxBN)xBDII 


NME 


NUT75 


(BDIIxBN)xBDII 


NME 


NUT110 


(BDIIxBN)xBDII 


NME 


NUT122 


(BDIIxBN)xBDII 


NME 


NUT129 


(BDIIxBN)xBDII 


NME 


NUT7 


(BDIIxSPRD)xBDII 


EAC 


NUT41 


(BDIIxSPRD)xBDII 


EAC 


NUT42 


(BDIIxSPRD)xBDII 


EAC 


NUT47 


(BDIIxSPRD)xBDII 


EAC 


NUT84 


(BDIIxSPRD)xBDII 


EAC 


NUT58 


(BDIIxSPRD)xBDII 


NME 


NUT68 


(BDIIxSPRD)xBDII 


NME 


NUT74 


(BDIIxSPRD)xBDII 


NME 


NUT89 


(BDIIxSPRD)xBDII 


NME 


NUT91 


(BDIIxSPRD)xBDII 


NME 



The tumour material was derived from cross between females from the EAC- 
susceptible BDII strain and males from BN or SPRD, respectively. NUT- tumour 
developed in the backcross (N1) progeny. 



cancer [16,17]. The first documentation of a miRNA 
abnormality in cancer was reported by Croce and col- 
leagues in 2002, who found that two miRNAs, miR-15 
and miR-16, were lost or down-regulated in most 
patients with chronic lymphocytic leukaemia (CLL) [18]. 
Furthermore, altered miRNA expression profiles have 
been reported in other cancers, such as colon [19], 
breast [20] and endometrial cancer [21]. Identification of 
deregulated miRNAs in different cancer types, and 
knowledge of their functional significance, will advance 




0 H 1 1 1 1 1 

1 2 3 4 5 6 

Number of genes 

Figure 1 Accumulated SD using Norm Finder, indicating five 

reference genes (black marker) as the optimal choice. 

\ ) 



our understanding of the aetiology of cancer and ultim- 
ately lead to development of valuable therapeutic tools. 

Quantitative polymerase chain reaction (qPCR) is one 
of the most important methods used for gene expression 
analysis, since it is highly specific and allows quantitative 
detection of small changes in gene expression. One of 
the most critical points in determination of a reliable 
expression pattern is removal of non-biological (experi- 
mental) variation from true biological variation. The fac- 
tors that influence the variability of qPCR expression 
values include the amount and quality of the starting 
material, the yields of the extraction process and the re- 
action efficiency. To deal with these factors it is import- 
ant to use a suitable normalization strategy in the data 
analysis. The commonly used option for normalization is 
to use stably expressed endogenous control genes [22]. 

A suitable endogenous control gene should have abun- 
dant and constant expression across all samples. I would 
be very advantageous to identify control genes that are 
stably expressed independent of cell type and tissue type. 
Since no such control gene has been identified, it is 
important to validate the expression of endogenous con- 
trol genes [13]. A good endogenous control for normali- 
zation of miRNA expression data must share similar 
properties with miRNAs, such as the size and stability of 
the RNA molecule. Some classes of small non-coding 
RNAs, not belonging to the miRNA class, are often 
expressed in an abundant and stable manner, making 
them good candidate control genes [14,15]. 

While many endogenous controls for miRNA expres- 
sion studies have been identified in different species 
such as human and pigs [23,24], no well-established en- 
dogenous controls for miRNA quantification have been 
identified in rat. One of the endogenous control genes 
that we identified as a suitable in this report, the U6 
snRNA, has been used as a reference gene in recently 
published studies on lung cancer [25]. 

In this study, we investigated the expression of 4.5S 
RNA (H) A, Yl, 4.5S RNA (H) B, snoRNA U87 and U6 
in 21 rat cell lines. NormFinder and geNorm were used 
to find suitable control genes from this set of candidates, 
and both algorithms identified the same five control 
genes as most stably expressed, only differing in their 
ranking order. This is not unexpected since the two 
algorithms use different mathematical models. GeNorm 
calculates the stability value defined as the arithmetic 
mean of all pair-wise comparisons, whereas NormFinder 
estimates both the overall expression variation as well as 
the variation between sample groups. Therefore, the re- 
sults of this study indicate that Yl should not be used as 
endogenous control since it is ranked as the least stable 
by both algorithms. It is noteworthy that the five 
most stable genes in this study are all expressed in the 
nucleus, while Yl is expressed in the cytoplasm. The 
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Student s £-test gives an additional clue to the stability 
of the endogenous control genes, where genes with the 
highest p values would be the most stable ones 
(Table 1). 

Beside the stability values, NormFinder also calculates 
the accumulated standard deviation, which determines 
the number of reference genes that should be used for 
normalization. The results from NormFinder indicate 
that use of the five most stable genes (U6, 4,5S RNA (H) 
A, U87, snoRNA and 4,5S RNA (H) B) provides the best 
normalization. However, using five reference genes 
would be both time consuming and expensive as well as 
require a lot of starting material. It may therefore be 
advantageous to use snoRNA and U87, which are highly 
ranked by both algorithms, as well as by the £-test. 

Conclusion 

In summary, careful selection of appropriate endogenous 
control genes plays a crucial role in expression studies. 
The candidate endogenous control genes should be vali- 
dated, because a gene that is stable in one cell line or tis- 
sue type is not necessarily stable in other cell lines/tissue 
types. We recommend use of U87 and snoRNA as con- 
trol genes in miRNA expression analysis in rat cells. 

Materials and methods 

Tumour material 

In a previous study females from the EAC-susceptible 
BDII strain were crossed to males from two other 
inbred rat strains that are not prone to develop EAC, 
BN/Han and SPRD-Cw3/Han (hereafter BN and SPRD). 
An Fl generation was generated and the Fl progenies 
were either backcrossed to BDII females to generate an 
Nl generation, or intercrossed to produce F2 progeny. 
Tumour formation was monitored by palpation. In 
cases of suspected tumour, the animals were sacrificed 
and a necropsy was done. Pathological analysis showed 
that the majority of the tumours that developed in the 
Nl progeny were classified as EAC, but in some cases 
no malignant cells were observed in the removed cell 
mass. These tissue samples were classified as non-/pre- 
malignant endometrium [26]. Spontaneously arising 
tumours developed in approximately 25% of the Fl, 
F2 and Nl offspring [27,28]. The NUT (backcross 



(Nl) uterine tumour)specimens represent non-/pre- 
malignant cells (NME) from the endometrium or 
endometrial adenocarcinomas (EAC) developed in the 
backcross (Nl) progeny. 

In the present study, a total of 20 NUT endometrial 
cell lines were studied, of which 10 were classified as 
NME and 10 as EAC. Ten of the cell lines were derived 
from the (BNxBDII)xBDII crosses and ten from the 
(SPRDxBDII)xBDII crosses (Table 2). A rat embryo 
fibroblast (REF) cell culture was used as normal control. 

Cell cultures 

The cell lines established from endometrial adenocarcin- 
oma and non-/pre-malignant tissues were cultured in 
Dulbecco s modified Eagle medium (DMEM) supplemented 
with 10% heat-inactivated fetal bovine serum, 100 IU/ 
100 ug ml - 1 penicillin/streptomycin, L-glutamine, MEM 
essential and non-essential amino acids, and MEM vitamins 
solution. The cells were grown at 37°C in an atmosphere of 
95% humidity and 5% C0 2 . The cells were harvested by 
trypsinization at 80-90% confluence (<1 x 10 6 cells). 

Candidate control genes 

A total of six candidate endogenous control genes were 
selected based on the availability of commercial primers. 
They are all small non-coding RNAs; five small nuclear 
RNAs (4.5S RNA (H) A, 4.5S RNA (H) B, snoRNA, U87 
and U6) and one small cytoplasmic RNA (Yl). Details 
regarding the selected RNAs are listed in Table 3. 

RNA extraction, reverse transcription and qPCR 

Total RNA including miRNA was isolated from the cell 
lines using a mirVana miRNA Isolation Kit (Ambion) 
following the manufacturer s protocol. Quality and quan- 
tity of the RNA samples were determined in a NanoDrop 
ND-1000 Spectrophotometer (NanoDrop Technologies, 
USA). All the RNA samples had a 260/280 absorbance 
ratio of 2.0-2.1. Total RNA samples were converted to 
cDNA and aliquoted into triplicates using TaqMan 
microRNA reverse transcription kit and TaqMan miRNA 
primers (Applied Biosystems). After the RT step, the real- 
time PCR reactions were performed according to the 
manufacturers instructions; the reaction consisted of 
1.33 (il cDNA (10 ng), 1.0 ul TaqMan Small RNA Assay 



Table 3 Candidate control genes selected for evaluation of expression stability 


Assay name 


Name 


Accession no 


Type of RNA 


Database 


4.5S RNA(H) A 


4.5S RNA(H) variant 5 


AY228151 


snRNA 


GenBank 


Y1 


Y1 scRNA gene 


U84683 


scRNA 


GenBank 


4.5S RNA(H) B 


4.5S RNA(H) variant 1 


AY228147 


snRNA 


GenBank 


snoRNA 


E2 small nucleolar RNA gene 


U64702 


snRNA 


GenBank 


U87 


small nucleolar RNA U87 


AF272707 


snRNA 


GenBank 


U6 


U6 small nuclear 1 (RNU6-1) 


NR_004394 


snRNA 


RefSeq 
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(20X), 10 |il TaqMan Universal PCR Master Mix II and 
7.67 ul nuc lease-free water. All reactions were performed in 
triplicates, including the no-template control (NTC). The 
reactions were run on an Applied Biosystems 7300 Real 
Time PCR system with the following thermal cycles: one 
cycle of 95°C for 10 minutes; 40 cycles with a denaturation 
step at 95 °C for 15 seconds and an annealing/extension step 
at 60°C for 60 seconds. 

Data analysis 

Relative quantities for each candidate gene were calcu- 
lated using the comparative Ct method. The Ct value 
represents the cycle number at which the fluorescence 
passes the defined threshold. One-way analysis of vari- 
ance (ANOVA) was used to test for any significant dif- 
ference among replicates (p < 0.05). Differences in gene 
expression between malignant and non-/pre-malignant 
samples were calculated by the Students £-test (p < 0.05). 

GenEx software (MultiD Analyses AB, Goteborg, 
Sweden) was used to analyze the stability of candidate 
genes with the geNorm and NormFinder algorithms. 
GeNorm calculates, for each candidate control gene, the 
pairwise variation with all other candidate genes as the 
standard deviation of the log-transformed expression ra- 
tios, followed by calculation of the M-value. The M-value 
describes the average pairwise variation for a particular 
control gene with all other candidate genes. At each com- 
parison the gene that shows largest variation to all the 
genes is eliminated. The process is repeated until there is 
only one pair of genes left. These last two genes are 
recommended as the optimum control genes [29]. 

NormFinder calculates the variance of gene expression 
both within groups and between groups (such as normal 
vs cancer). The intra- and intergroup variances are com- 
bined into a single stability value and the algorithm 
ranks the candidate control genes according to this com- 
bined variance. The main differences between geNorm 
and NormFinder are that NormFinder takes into 
account both the inter- and intragroup variation [30]. 

Abbreviations 

miRNA: microRNA; EAC: Endometrioid adenocarcinoma; NME: NonVpre- 
malignant cells; qPCR: Quantitative polymerase chain reaction. 
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